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Magnetization Plateaus in One Dimensional S = 1/2 Heisenberg Model 
with Dimerization and Quadrumerization 
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The one dimensional 5 = 1/2 Heisenberg model with dimerization (l—j) and quadrumerization 
(5) in the magnetic field is studied by means of the numerical exact diagonalization of finite size 
systems and the conformal field theory. It is found that the magnetization plateau at half of 
the saturation value exists for 5^0. For (5 = 0, this model is described by the conformal 
field theory with central charge c = 1 at this value of magnetization. The critical exponent u 
which characterizes the ^-dependence of the width of the plateau is calculated using the level 
spectroscopy method. The ji-dependence of the critical exponent f is found to be non-monotonic 
and discontinuous at j =0. The effective theory of the magnetization plateau is also presented 
for various limiting cases. 

KEYWORDS: Heisenberg chain, spin gap, exact diagonalization, magnetization plateau, conformal field theory, 
central charge, scaling dimension, dimer gas model. 



§1. Introduction 

Recently, the Heisenberg chains with a variety of spa- 
tial structures have attracted great interest. Although 
the uniform S = 1/2 ajitiferromagnetic Heisenberg chain 
can be solved exactlyEP and its ground state is known as 
a gapless spin liquid, various spatial structures such as 
spin-Peierls dimerization, zig-zag chain or ladder struc- 
ture can drive this ground state into spin gap states. 

One the other hand, the magnetization plateaus in one 
dimensional Heisenberg. |ChAi|nsJpye-.bejS^ 
by a lot of researcherffliyBy-ElByy yBy y y) quite 
recently. This state can be regarded as the field induced 
spin gap state in which spins are partly quenched by 
the magnetic field and remaining spins form a state with 
a finite energy gap i.e. the magnetization plateau. It is 
evident that such magnetization plateaus can be realized 
in various non-homogenecuis Heisenberg chains. 

In the previous work,E3' two of the present authors 
(Chen and Hida) investigated the ground state phase 
diagram of the S = 1/2 Heisenberg model with coex- 
isting dimerization and quadrumerization. Similar poly- 
merized Hpitynj-iprcr rhajng are also Studied by many au- 
thorsE3'oE3E3'U'E3E3'E3'' In this paper, we further in- 
vestigate the magnetization plateaus in this model using 
the numerical exact diagonalization method and the con- 
formal field theory. 

This paper is organized as follows. In the next sec- 
tion, the model Hamiltonian is defined. The numerical 
results for the magnetization plateaus are presented in 
§3. Analyzing the exact diagonalization data by the con- 
formal field theory, we calculate the central charge c and 
the critical exponent v characterizing the opening of the 
plateau. The effective theory in several limiting cases is 



also pesented. The last section is devoted to summary 
and discussion. 

§2. Model Hamiltonian 

The Hamiltonian of the one dimensional dimerized and 
quadrumerized = 1/2 Heisenberg chain in the mag- 
netic field is given by 



2N 



2N 



n = jJ2 S21-1S21 + ^(1 + i-iy-'s)S2iS2i+i 

1=1 1=1 

4N 

-gii^HY^St, (2.1) 



1=1 



where 1 — j(— 00 < j < 00) and i5(— 1 < S < 1) rep- 
resent the degree of dimerization and quadrumerization, 
respectively. The magnetic field, the electronic g-factor 
and the Bohr magneton are denoted by H, g and /ib , re- 
spectively. In the following, we take the unit g/xe = 1 and 
5 > without loss of generality and assume the periodic 
boundary condition(5^4Ar+i = S'l). For j > 0, this model 
can be regarded as the quadrumerized antiferromagnetic- 
antiferromagnetic alternating chain while for j < 0, it 
can be regarded as the quadrumerized ferromagnetic- 
antiferromagnetic alternating chain. This model ap- 
proaches the dimerized spin-1 antiferromagnetic Heisen- 
berg chain in the limit j — 00. According to the 
criterion by Oshikawa et al.,EP it is possible that the 
present system exhibits the magnetization plateau at 
= 1/4 where is the magnetization per site de- 



fined by 77?/ 



1 

4Ar 



AN 

1=1 



§3. Numerical Results 
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3.1 Magnetization Plateaus 

We use the exact diagonaUzation method to calculate 
the magnetization curve of the present model. Typically 
the magnetization curve has the form shown schemat- 
ically in Fig. 1 with a plateau at = 1/4. The 5- 
dependence of the 4 critical fields HcQ, Hcli Hc2 and Hs 
are shown in Fig. 2 for (a)j = 5, (b) 0.5, (c) —0.5 and (d) 
—5. The critical field iJ^o corresponds to the energy gap 
in the absence of magnetic field, Hqi to the energy gap 
between the lowest energy with — N and that with 
— N —1, Hc2 the energy gap between the lowest en- 
ergy with M"" = N + 1 and that with M"" = N and is 
the saturation field. Here is the total magnetization 
defined by = ANm^. The magnetization plateaus 
at = 1/4 open between the Hd and Hc2 as shown 
in the Fig. 2(a-d) by the symbols o and o, respectively. 
The values |af Hd and iJc2 in Fig- 2 are obtained by 
the Shanks'u-' extrapolation of the results of the Lanc- 
zos exact diagonalization for 4N — 12, 16, 20, 24 and 28. 
The points at which iJ^o vanishes are determined using 
the method of Ref. |l^. The symbol • represents H^o 
obtained similarly from the data for 4N — 8, 12, 16, 20 
and 24. The thick solid line is the saturation field Hs 
which can be calculated analytically as follows. 

The saturation field is given by the energy difference 
AE between the loweKLenergy with = 2N — 1 and 
that with Af^ = 2A^.dc3-' The latter is simply given by 
the energy of the fully polarized state. The former is the 
plane wave state of a single inverted spin. The energy 
difference AE is given by the lowest eigenvalue of the 
matrix, 



l-<5 
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l-i5 „i4fc 
2 ^ 



1+5 
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2 



1-S, 



-iik 



1+S 
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1+S 
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1-S 
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(3.1) 

Where the wave number is denoted by k. We find that 
AE is the lowest for fc = 0. Thus, the saturation field 
Hs is obtained as. 



Hs 



-AE 



ij + 2) + + 4^2 



(3.2) 



In view of these figures, it is quite likely that the 
plateau at = 1/4 opens for infinitesimal 5. There- 
fore we expect that 5 = is the critical point. It should 
be also noted that the similar plateau is found in the 
S = 1 dimerized antiferromagnetic Heisenberg chain not 
only theoreticallyO'EP but also experimentally .Er This can 
be regarded as the j ^ ~oo limit of the present model. 
Comparing Fig. 2(a,b) and (c,d), the width of the mag- 
netization plateaus glows more rapidly for j > Q than 
for j < 0. This suggests that the critical exponent v of 
the plateau width defined by _ffc2 — H^i ~ 5^ behaves 
differently depending on the sign of j. In the next sub- 
section, we g|tu.dy this point in more detail by the level 
spectroscopyE3 analysis of the Hamiltonian with 5 — 
based on the conformal field theory. 



3. 2 Critical Exponent v of the Magnetization Plateau 

Assuming that the ground state with = 1/4 is crit- 
ical and conformally invariant at (5 = , we estimate the 
central charge c from the Lanczos exact diagonalization 
data. It is known that the finite size correction to the 
ground st at^^ ep ergy is related with the central charge as 
follows, 



J-EJN,M^ 
2N ' 



TT 



1 



6 '(2iV)2 



(3.3) 



where Vs is the spin wave velocity, Eg{N, M^) is the 
ground state energy with system size AN and magne- 
tization and e(m^) is the ground state energy per 
unit cell in the thermodynamic limit with corresponding 
value of . It should be noted that the size of the unit 
cell is 2 for (5 = 0. In order to determine Vs, we also 
calculate the excitation energy with fixed wave vector k 
under the periodic boundary condition with = N. 
The system sizes are AN = 12, 16, 20, 24 and 28. We de- 
note the lowest energy with fixed k by Ek{N, M^) and 
the lowest one among them is the ground state energy 
Eg{N, M^). For j > 0, the ground state has fc = ir- 
respective of the value of M^. For j < 0, on the other 
hand, the grpund state has fc = for even and fc = tt 
for odd M^@) 

Thus, the velocity Vs is estimated by 



2N 

= [E,AN,An 

2,71 



Eg{N,M^)], (3.4) 



where fci is equal to 2tt/{2N) for j > and j < with 
even M^, whereas it is tt — 27r/ (2iV) for j < with odd 

From Eq. (3.3), the gradient A of the plot 
£;g(iV, M^)/2N versus l/(2iV)2 is equal to 7rcWs/6 which 
has negligible size dependence. The finite size central 
charge c is determined as 6^/(7rws). The extrapola- 
tion procedure for c is shown in Fig. 3 which yields 
c ~ 1.019 ± 0.001 for j = 0.5. Figure 4 shows the j- 
dependence of the central charge c which indicates that 
c is close to unity everywhere. The error bars due to 
the extrapolation procedure are less than the size of the 
symbols. Thus we may safely assume that the critical 
line 5 = with — 1/4 is described by the Gaussian 
model with conformal charge c = 1 given by, 



«G = i/<i.[...AW + |(|rt. 



(3.5) 



where (j) is the boson operator compactified as < < 
\/2i: and 11 is the momentum density conjugate to (j) 
which satisfies [(/)(x), n(x')] = i5{x — x'). 

The scaling dimension Xn of an operator On is related 
with the energy eigenvalue En{N, M^) of the wta.t^p gpy 
ated by applying the operator to the ground state 

0iS • X Ti 



lim Xn{N), where 

N—^oo 



Xn{N) 



2N 
2nvs 



[E„iN,M^)-EgiN,M')]. (3.6) 



We can identify the correspondence between the opera- 
tors in boson representation and the eigenstates of the 
spin chains by comparing their symmetry O'EZP In this 
way, from the excitation energies of the eigenstates cor- 
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responding to the operators Oi = cos\/20 and O2 — 
sin ^/2(j), we can determine their scaling dimensions xi 
and X2 using the relation (3.6). On the other hand, 
both xi and X2 should be equal to K/2 as determined 
from their correlation function. Thus the value of K is 
determined from the numerically obtained values of xi 
and/or X2. Actually, to reduce the finite size correction 
to 0{1/N'^), it is more convenient to use the combina- 
tion, 



K{N) = xi{N) + X2{N), 

as proposed by Kitazawa and NomuraH-* 
the formula 



K{N) + 



ci 



C2 



(27V)2 



(3.7) 

We assume 

(3.8) 



(27V)4 

to extrapolate K{N) to — > oo. The extrapolation 
procedure is shown in the Fig. 5 which gives K = 
1.5227 ± 0.0003 for j = 0.5. The error is estimated from 
the difference between the values extrapolated from the 
data for AN = 12, 16, 20, 24 and AN = 16, 20, 24, 28. 

From the symmetry consideration, the effect of 
quadrumerization can be represented by the term pro- 
portional to cos v^i/) which explicitly breaks the under- 
lying U{1) symmetry. Such a term can be also explicitly 
derived iDiLjpi? = 1 /4 by the usual Abelian bosonization 
scheme .E3E3'I3'I3eP Thus the low energy properties of our 
model can be described by the sine-Gordon model given 

by, 



2tt 
27ra2 



dx 



K^dx' 



dx cos 



(3.9) 



where the a is the lattice constant. The second term 
of the Hamiltonian generates the energy gap. 

From the renormalization group theory and conformal 
field theory, the critical exponent of the energy gap v 
produced by the cos-\/20 term is related to its scaling 
dimension xi as 



1 



1 



2 -XI 



K 
2 



(3.10) 



Figure 6 shows the variation of the critical exponent v 
with j. The error bars due to the extrapolation proce- 
dure are less than the size of the symbols. For j > 0, 
the critical exponent v is close to 0.8 and slightly de- 
pends on j with a weak maximum i/ ~ 0.81 at j = 1, 
while for j < 0, the critical exponent v increases with j 
and appraoches 2 as j —^ —0. The positive values of v 
ensure that the point (5 = is the critical point (unsta- 
ble fixed point) and the plateau opens for infinitesimal 
S. These values are consistent with the correlation func- 
tion exponents ri{— 1/|£C) and ri^{— K) for j = —2 and 
0.5 obtained by SakaLQ The values for large negative 
j are also consistent with the values of 7] and f]^ for the 
S — \ antiferromagn£±ic Heisenberg chain obtained by 
Sakai and TakahashiH which give v = 1.74 ± 0.1. Ob- 
viously, V — 1 ior j — Q. It is remarkable that there is a 
discontinuity in at j = 0. 



3. 3 Analytic Results for the Critical Exponent v 
3.3.1 Undimerized case : j = 1 



For j = 1 and 6 



the solution by the Bethe's 



hypothesis is available. a In this case, the parameter 
K is determined from the dressed charge to yield v ~ 
0.81 for = 1/4 by Fledderjohann et al!3i3 which is 
consistent with our numerical results. 
3.3.2 Mapping onto the effective unmagnetized XXZ 
chain for j ~ and j » 1 

The magnetized state of the present model on the 
plateau can be mapped onto the unmagnetized state of 
the effective XXZ chain in the limiting cases j ^ and 
j ~* oo using the method proposed by Totsuka.EP 

For small values of j, the dimers are formed on the 
1 + ^-bonds and 1 — J-bonds in the ground state, if the 
magnetic field is absent. When the magnetic field is 
switched on, some of the dimerized bonds are broken 
and change into the |tT> state. Other states on the 
1 ± ^-bonds are energetically unfavourable and can be 
discarded. Then the present system can be described as 
a gas of the interacting dimers as follows 



T-i' = j Y^[n^n,+i - {dfd,+i + d+^^di)] 



2N 



1=1 
2^ 



2JV 



^(if-l)n, + <5^(-l)*n„ (3.11) 



4=1 



4=1 



where df and di are the fermion operators and = 
dfdi- The state with ni=l corresponds to the state in 
which the z-th 1 ± (5-bond is occupied by a dimer, while 
ni=0 to the |tT> state. This Hamiltonian can be further 
mapped onto the 5' = 1/2 XXZ chains by the Jordan- 
Wigner transformation with the correspondence 

i (3.12) 



— S/ 



where Si is the effective spin-1/2 operator, 
effective XXZ Hamiltonian is given by 



Then, the 



cfT 



2N 

■E 

i=l 



2N 

E( 

1=1 



H - l)S^ 



2N 
i=l 



where A 



(3.13) 

1/2 for j > and A = -1/2 for j < 0. It 
should be noted that the staggered field term propor- 
tional to 6 appears because the energy required to break 
the dimer bond is different according whether i is even 
or odd. 

The parjEmfiter K can be calculated by the Bethe's 
hypothesisES'El^ as 

Therefore we have K = 3/2 for j ^ -|-0 and AT = 3 
for j —0. In the bosonized language, the staggered 
field term corresponds to sin_\/20 term.c3 Therefore we 
can again use the formula ( |3.10| ) to obtain v is 4/5 for 
j +0 and 1/ = 2 for J ^ —0. 

On the other hand, for j — > +oo, the dimers reside on 
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the j-bonds in the ground state. Using the similar argu- 
ment as the small j case, the effective XXZ Hamiltonian 
can be written as 

Ks - Ei 2 ^^^^^^^ + ^'^'+1 + 

i=l 



2N 



(3.15) 



with A = 1/2 resulting in K — 3/2. In this case, the 
dimerization term appears in the effective Hamiltonian 
which has the bosonized form cos v^^Bj-* This term also 
has the same scaling dimension K /2. Therefore we have 
again = 4/5 for j — > oo. Thus the analytic results 
explain our numerical results quite well. 

§4. Summary and Dicussion 

The magnetization plateaus at — 1/4 in the one di- 
mensional S* = 1/2 Heisenberg model with dimerization 
and quadrumerization are investigated by the exact diag- 
onalization of the finite size systems. The central charge 
at the critical point (5 = is calculated to be close to 
unity. Under the assumption of the universality class of 
the c = 1 Gaussian model, thus, the critical exponent 
of the plateau width v is obtained by the level spec- 
troscopy method from the numerically calculated energy 
spectrum. The critical exponent v is almost constant 
and close to 0.8 for j > 0, while for j < 0, it strongly 
depends on j and increases up to 2 as j approaches —0 
as shown in Fig. 6. 

For j — > ±0 and j +oo, the present model is 
mapped onto the effective unmagnetized XXZ model 
based on the dimer gas model. This gives the results 
consistent with the numerically obtained values of v. At 
j ~ I, the result is also consistent with the analytical 
calculation using the Bethe's hypothesis. 

According to the present calculation, the magnetiza- 
tion plateau at — 1/4 should be always observed 
in the presence of quadrumerization. Therefore it must 
be interesting to synthesize the linear chain material 
which has the present structure. It should be also re- 
marked that it is easier to observe the plateau for posi- 
tive j rather than negative j as far as the strength of the 
quadrumerization is weak. 

The value of ly is also related to the exponent of the 
energy gain AE due to quadrumerization at small S as 
AE (X 5'^". For j > 0, 2j/ - 1.6 < 2. Therefore the 
energy gain due to quadrumerization is larger than the 
lattice deformation energy and quadrumerization takes 
place spontanuously at — 1/4 in the similar way as 
the spin-Peierls instability. Actually, such instability can 
take place for arbitrary values of magnetization as far 
as v < 1. In the case of quadrumerization, the wave 
number of the spatial modulation is commensurate with 
the original lattice. Therefore the lattice distortion is 
pinned to the lattice and the plateau will be stabilized. 
Observation of such field induced spin-Peierls transition 
would be also an interesting problem from the experi- 
mental side. On the other hand, for j < 0, v is always 
larger than unity and no spontanuous quadrumerization 



is expected. 
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Fig. 1. Schematic magnetization curve of the present model. 



Fig. 2. The critical fileds for (a)j = 5, {h)j = 0.5, {c)j = -0.5 
and (d)jr' = —5. The critical fields J^qQ' ^cl ^c2 '='"^ "^P" 
resented by •, o and o, respectively. The solid thick line is the 
saturation field i?s- The thin solid lines are just guides for eye. 



Fig. 3. The extrapolation procedure of the finite size central 
charge c for j = 0.5. 



Fig. 4. The j dependence of the numerically obtained central 
charge c. 



Fig. 5. The extrapolation procedure of K at j = 0.5. 



Fig. 6. The j dependence of the critical exponent v. 
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Fig. 2(a) Chen, Hida, Nakano 




Fig. 2(b) Chen, Hida, Nakano 




Fig. 2(c) Chen, Hida, Nakano 




Fig. 2(d) Chen, Hida, Nakano 
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